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c/5 ; Abstract 

*r^ • We present a detailed analytical and numerical study for the spreading 

'^ ■ of infections in complex population networks with acquired immunity. We 

show that the large connectivity fluctuations usually found in these networks 
O I strengthen considerably the incidence of epidemic outbreaks. Scale- free net- 

works, which are characterized by diverging connectivity fluctuations, exhibit 
\^J ■ the lack of an epidemic threshold and always show a finite fraction of infected 

f-~^ ■ individuals. This particular weakness, observed also in models without im- 

^^ . munity, defines a new epidemiological framework characterized by a highly 

f->^ i heterogeneous response of the system to the introduction of infected individ- 

^D ' uals with different connectivity. The understanding of epidemics in complex 

networks might deliver new insights in the spread of information and diseases 
in biological and technological networks that often appear to be characterized 
by complex heterogeneous architectures. 
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I. INTRODUCTION 

The epidemiology of heterogeneous networks has largely benefitted from the need of 
understanding the spreading of human sexual diseases in the complex web of sexual part- 
nership [|l]-0]. Epidemic modeling considered that population groups can be characterized 
in classes having different sexual activity or number of sexual contacts. This fact leads to 
models dealing with heterogeneous populations which are known to enhance the spread of 
infections as well as make them harder to eradicate (for a review see [^). In this perspec- 
tive, a limiting case is represented by the newly identified classes of complex networks (for 
a review see ) . The highly heterogeneous topology of these networks is mainly reflected 
in the small average path lengths among any two nodes (small- world property) |P,0, and in 
a power law distribution (scale-free property), P{k) ~ k~'^~"' , for the probability that any 
node has k connections to other nodes !§]. While regular networks present finite connectiv- 
ity fluctuations {{k) ~ (fc^)), scale-free (SF) networks are a limiting case of heterogeneity 
where connectivity fluctuations are diverging if < 7 < 1. In other words, the network 
nodes possess a statistically significant probability of having a virtually unbounded number 
of connections compared to the average value. SF networks find real examples in several 
technological networks such as the Internet P,|TO| and the world- wide- web (WWW) [0], as 



well as in natural systems such as food- webs, and metabolic or protein networks [^]. The 
need to understand the dynamics of information transmission, the error tolerance [p!^-[T^ 
and other properties of complex networks has therefore called for the study of epidemic 
modeling in complex networks |[T5|-|T9|. 

A surprising result, originated by the inspection of the susceptible-infected-susceptible 
(SIS) model, has shown that the spread of infections is tremendously strengthened on SF 
networks [|18|,|19|. Opposite to standard models, epidemic processes in these networks do not 
possess any epidemic threshold below which the infection cannot produce a major epidemic 
outbreak or an endemic state. In principle, SF networks are prone to the persistence of 
diseases whatever infective rate they may have. This feature reverberates also in the choice 



of immunization strategies ||20|-p2| and changes radically many standard conclusions on 
epidemic spreading. This study appears particularly relevant in the case of technological 
networks, for instance for the spreading of digital viruses in the Internet [1^, and it has soon 
been generalized by showing that also the susceptible-infected-removed (SIR) model show 
the same absence of epidemic threshold ||2^. These results highlight the study of epidemic 
models in complex networks as potentially relevant also in human and animal epidemiology 
p3| , as confirmed recently by the experimental observation that the web architecture of 
sexual contacts is best described by a scale-free topology in which individuals have widely 
different connectivities [P^]. 

In this paper we provide a detailed analytical and numerical study of the SIR model 
on two prototype complex networks: the Watts-Strogatz (WS) model and the Barabasi - 
Albert (BA) model. The first model is a small-world network with bounded connectivity 
fluctuations, while the second one is the prototype example of SF network. The analytical 
approach allows us to recover the total size of the epidemics in an infinite population, in 



agreement with earlier estimates [2^. We are able to find the analytic expression for the 
critical threshold as a function of the moments of the connectivity distribution and we 
confirm the absence of any finite threshold for connectivity distributions P{k) ~ /c~^~^ with 



< 7 ^ 3. We obtain the general analytic expression for the total density of infected 
individuals and the epidemic threshold at arbitrary 7 values. Finite size network effects 
can be easily evaluated from the analytic expressions. Time evolution and other effects of 
heterogeneity such as the relative infection incidence in different connectivity classes can 
be predicted. In order to confirm the analytical findings we perform large scale numerical 
simulations on the WS and BA networks. Numerical results are in perfect agreement with 
the analytical predictions and confirm that the interplay of complex networks topology and 
epidemic modeling leads to a new and interesting theoretical framework, whose predictions 
and implications need to be exhaustively explored. 

During the completion of this paper we became aware of a preprint by Lloyd and May 
| 25| which reports a comprehensive study of the SIR model in scale free networks. This work 



extends the preliminary account provided in Ref. EH 



II. THE SIR MODEL 

Our theoretical understanding of epidemic spreading is based on compartmental models, 
in which the individuals in the population are divided in a discrete set of states @,|^ . In this 
framework, diseases which result in the immunization or death of infected individuals can 
be characterized by the classical susceptible-infected-removed (SIR) model @J2^. In this 
model individuals can only exist in three different states: susceptible (healthy), infected, or 
removed (immunized or dead). In a homogeneous system, the SIR model can be described 
in terms of the densities of susceptible, infected, and removed individuals, S{t), p{t), and 
R{t), respectively, as a function of time. These three magnitudes are linked through the 
normalization condition 

S{t)+p{t) + R{t) = l, (1) 

and they obey the following system of differential equations: 

^ = -/ip + XkpS, (2) 

dR 

11 = ^P- 

These equations can be interpreted as follows: infected individuals decay into the removed 
class at a rate p, while susceptibles individual become infected at a rate proportional to both 
the densities of infected and susceptible individuals. Here, A is the microscopic spreading 
(infection) rate, and k is the number of contacts per unit time that is supposed to be constant 
for the whole population. In writing this last term of the equations we are assuming the 
homogeneous mixing hypothesis |Q, which asserts that the force of the infection (the per 
capita rate of acquisition of the disease by the susceptible individuals) is proportional to the 
density of infectious individuals. The homogeneous mixing hypothesis is indeed equivalent 
to a mean-field treatment of the model, in which one assumes that the rate of contacts 
between infectious and susceptibles is constant, and independent of any possible source of 



heterogeneity present in the system. Another imphcit assumption of this model is that the 

time scale of the disease is much smaller than the lifespan of individuals; therefore we do 

not include in the equations terms accounting for the birth or natural death of individuals. 

The most significant prediction of this model is the presence of a nonzero epidemic 



threshold Ac l^g. If the value of A is above Ac, A > Ac, the disease spreads and infects a 
finite fraction of the population. On the other hand, when A is below the threshold, A < Ac, 
the total number of infected individuals (the epidemic incidence), -Roo = linij^oo -R(^), is 
infinitesimally small in the limit of very large populations (the so-called thermodynamic 
limit ||2^). In order to see this point, let us consider the set of equations (^, in which, 
without lack of generality, we set /x = 1. Integrating the equation for S{t) with the initial 
conditions -R(O) = and 5'(0) ^ 1 (i.e., assuming p(0) ~ 0, a very small initial concentration 
of infected individuals), we obtain 

^(t)=e-^^^W. (3) 

Combining this result with the normalization condition (|I|), we observe that the total number 
of infected individuals -Rqo fulfills the following self-consistent equation: 

i?oo = l-e^^^^°°. (4) 

While Roo = is always a solution of this equation, in order to have a nonzero solution the 
following condition must be fulfilled: 



dR_ 

This condition is equivalent to the constraint A > Ac, where the epidemic threshold Ac takes 
the value Ac = A; ~^ in this particular case. By using a Taylor expansion at A ~ Ac it is 
then possible to obtain the epidemic incidence behavior Roo ~ (A — Ac) (valid above the 
epidemic threshold). It is worth remarking that the threshold mechanism is related to the 
basic reproductive rate TZq ~ Xk (not to be confused with removed individuals) usually 
considered by epidemiologist. Only if TZq is larger than unity the infection can sustain itself, 



obviously defining a threshold in the spreading rate A |^ . As well, in the language of the 



physics of nonequilibrium phase transition p^, the epidemic threshold can be considered 



as completely equivalent to a critical point. In analogy with critical phenomena, we can 
consider i?oo as the order parameter of a phase transition and A as the tuning parameter. 
In particular, it is easy to recognize that the SIR model is a generalization of the dynamical 
percolation model, that has been extensively studied in the context of absorbing-state phase 
transitions ||2^ . 



III. THE SIR MODEL IN COMPLEX NETWORKS 

In order to address the effects of contact heterogeneity in epidemic spreading, let us 
consider the SIR model defined on a network with general connectivity distribution P{k) 
and a finite average connectivity (k) = Y^k^^i^)- Each node of the network represents an 
individual in its corresponding state (susceptible, infected, or removed), and each link is a 



connection along which the infection can spread. The disease transmission on the network 
is described in an effective way: At each time step, each susceptible node is infected with 
probability A, if it is connected to one or more infected nodes. At the same time, each 
infected individual becomes removed with probability /i, that, without lack of generality, we 
set equal to unity. 

In order to take into account the heterogeneity induced by the presence of nodes with 
different connectivity, we consider the time evolution of the magnitudes pfc(t). Skit), and 
Rkit), which are the density of infected, susceptible, and removed nodes of connectivity k at 
time t, respectively. These variables are connected by means of the normalization condition 

Pk{t) + Skit) + Rkit) = l. (6) 

Global quantities such as the epidemic incidence are therefore expressed by an average over 
the various connectivity classes; i.e. i?(t) = I]fc-P(^)-Rfc(^)- ^^ the mean-field level, these 
densities satisfy the following set of coupled differential equations: 



dpkit) 

dt 
dSkit) 

dt 
dRkit) 



-Pkit) + \kSkit)eit), (7) 

-XkSkit)Qit), (8) 

Pkit). (9) 



dt 
The factor G(t) represents the probability that any given link points to an infected site. This 



quantity can be computed in a self-consistent way [|TB[: The probability that a link points to 
a node with s links is proportional to sPis). Thus, the probability that a randomly chosen 
link points to an infected node is given by 

ro.u^ _ ^kkPik)pkit) _ j:kkPik)pkit) .._. 

^^" EsSPis) - {k) ■ ^'"^ 

In this approximation we are neglecting the connectivity correlations in the network, i.e., 
the probability that a link points to an infected node is considered independent of the 
connectivity of the node from which the link is emanating. A more refined approximation 
would consider the network correlations as given by the conditional probability Pik/k') 
that a node with given connectivity k' is connected to a node with connectivity k p9| . 



Nevertheless, as we will see in the next sections, this rather crude approximation is quite 
able to give account of many of the properties shown by computer simulations of the model. 
The equations (0)-®, combined with the initial conditions -Rfc(O) = 0, Pfc(O) = p°, and 
'S'fc(O) = 1 — p° completely define the SIR model on any complex network with connectivity 
distribution Pik). We will consider in particular the case of a homogeneous initial distribu- 
tion of infected nodes, p° = p°. In this case, in the limit p° — *> 0, we can substitute pfc(O) — 
and S'fc(O) ~ 1. Under this approximation, Eq. (|]) can be directly integrated, yielding 

Skit) = e-^^'^W (11) 

where we have defined the auxiliary function 



(/)(t) = /* e{t')dt' = 7^ E kP{k)R,{t), (12) 

and in the last equality we have made use of the definition (|TD|). It is worth remarking 
that the above equations are similar to those obtained in the case of HIV dynamics in 
heterogeneous populations 

In order to get a closed relation for the total density of infected individuals, it results more 
convenient to focus on the time evolution of the averaged magnitude </>. To this purpose, let 
us compute its time derivative: 

^ = ^ E kP{k)p,{t) = ^ E kP{k){l - R,{t) - S,{t)) (13) 

= l-^(t)-^J2kP{k)Sk{t). (14) 

\'^/ k 

Introducing the obtained time dependence of Sk(t) we are led to the differential equation 
for 0(t) 

^ = l-^it)-^Y.kP{k)e-'''''^. (15) 

Once solved Eq. ([T5|), we can obtain the total epidemic incidence i?oo as a function of 
0OO = linit-^oo <P(t)- Since Rk{oo) = 1 — Sk{oo), we have 

i?o, = E^(^)(l-e-"'^°°)- (16) 

k 



Equations (|T5|) and ( [T6| ) constitute thus an alternative representation of the model, with 
respect to Eqs. (D-(D. 

For a general P{k) distribution, Eq. (|TB|) cannot be solved in a closed form. However, we 
can still get useful information on the infinite time limit; i.e. at the end of the epidemics. 
Since we have that Pk{oo) = and consequently limt^oo d4'{t) / dt = 0, we obtain from 
Eq. (|15D the following self-consistent equation for 0oo 

'^oo = l-7^E^^We"''^°°- (17) 

\'^/ k 

The value 0oo = is always a solution. In order to have a non-zero solution, the condition 

d 



(i-(i)i:*me— ») 



> 1 (18) 

<j)oo=0 



must be fulfilled. This relation implies 



-^^j:kpmxk) = x^-^>i. (19) 



This condition defines the epidemic threshold 



below which the epidemic incidence is null, and above which it attains a finite value. That 
is, the threshold is inversely proportional to the connectivity fiuctuations (fc^). For regular 
networks, in which (fc^) < oo, the threshold has a finite value and we are in the presence 
of a standard phase transition. On the other hand, networks with strongly fiuctuating 
connectivity distribution, show a vanishing epidemic threshold for increasing network sizes; 
i.e. (k"^) — > oo for A^ — > oo. The absence of any intrinsic epidemic threshold in this network 
can be understood by noticing that in heterogeneous systems the basic reproductive number 
7^0 contains a correction term linearly dependent on the fiuctuations (standard deviation) of 
the nodes' connectivity distribution [^p3[ . In SF networks the divergence of the connectivity 
fiuctuations leads to an 7?.o that always exceeds unity whatever the spreading rate A is. This 
ensures that epidemics always have a finite probability to survive indefinitely. It is worth 
remarking that real networks have always a finite size N and thus an effective threshold, 
depending on the magnitude of (k) and (/c^), that can be easily calculated as a function of 
A^. This apparent threshold, however is not an intrinsic quantity and it is extremely small 
for systems with large enough N. 



IV. EXPONENTIALLY DISTRIBUTED NETWORKS: THE WATTS-STROGATZ 

MODEL 

The class of exponential networks refers to random graph models which produce a con- 
nectivity distribution P{k) peaked at an average value {k) and decaying exponentially fast 
for k ^ (k) and k <^ (k). Typical examples of such a network are the random graph 
model ^^ and the small- world model of Watts and Strogatz (WS) [^. The latter has 
recently been the object of several studies as a good candidate for the modeling of many 
realistic situations in the context of social and natural networks. In particular, the WS 
model shows the "small- world" property common in random graphs [^]; i.e., the diameter 
of the graph — the shortest chain of links connecting any two vertices — increases very slowly, 
in general logarithmically with the network size [^. On the other hand, the WS model 



has also a local structure (clustering property) that is not found in random graphs with 



finite connectivity ||3T1 , |32[| . The WS graph is defined as follows ||3T|j3^: The starting point 
is a ring with N nodes, in which each node is symmetrically connected with its 2m nearest 
neighbors. Then, for every node each link connected to a clockwise neighbor is rewired to 
a randomly chosen node with probability p, and preserved with probability 1 — p. This 
procedure generates a random graph with a connectivity distributed exponentially for large 
k |^J5^, and an average connectivity (k) = 2m. The graph has small-world properties and 
a non-trivial "clustering coefficient"; i.e., neighboring nodes share many common neighbors 



3I|,|52I]. The richness of this model has stimulated an intense activity aimed at understand- 



ing the network's properties upon changing p and the network size N [|6|, p!5| , pl| -|3^ . At the 



same time, the behavior of physical models on WS graphs has been investigated, including 
epidemiological percolation models pr3|JT5|JT6|] and models with epidemic cycles [^ . 



In the following we focus on the WS model with p = 1; it is worth noticing that even 
in this extreme case the network retains some memory of the generating procedure. The 



network, in fact, is not locally equivalent to a random graph in that each node has at least K 
neighbors. In the limit p — *> 1, the connectivity distribution of the WS network, as defined 



previously, takes the form ||3^ 

P{k) 



m 



k—m 



for k > m. 



(21) 



{k — m)\ 

This is a Poisson distribution, with finite moments. Defining the factorial moments [^] 

(XO^ = (X(X - 1)(X - 2) ■ ■ ■ (X - r + 1)) , (22) 

we have for the distribution (pi|) 

{{k - mY) J = m"^ . (23) 

, the first moments of the connectivity distribution are given by 

(24) 
(25) 

(26) 



In particular, from Eq. 



(k) = 2m, 

kA = m(l + 4m). 



t 



and, in general. 



m(l + 6m + 8m 



(r) ~ (2m)' 



(27) 



for large m. 

For general regular networks, for which (A;") < oo for all values of n, Eqs. ([T5| ) and (|16D 
can be approximately solved in the limit 0(t) -^ 0, by expanding the exponentials under the 
summation signs. Thus, for the case of the total epidemic incidence i?oo in Eq- (P^, 



^oo ^^Pik)Xk(j)^ = {k) A0OO- 

k 



(28) 



That is, Roo is linearly proportional to 0oo- 

On the other hand, by expanding the exponential in Eq. (|TB|) and keeping the most 
relevant terms, we yield : 



d4> 

-- ~ 1 
dt 



Y.kkP{k){l - \k(t) + \^k^(p^ /2 



■1 + A 



(k^) ,.^jk^) 



(k) ^V ^ ' "(k) "^2{k))- 

The resulting previous equation can be exactly solved, to yield 

2(A - Ae) 1 



(29) 



m 



(P) (A;3) A2 + Ae-(^-^-)*A= 



(30) 



where Ac is defined as in Eq. (pO|). That is, for A < Ac, 0oo — ^ 0, while from A > Ac, we 
recover the well-known mean-field behavior (p^x, ~ (A — Ac), which translated to the total 
epidemic incidence i?oo yields 



8 



i?oo ~ (A - Ae) (31) 

In the particular case of the WS networks, we expect to observe the behavior dictated 
by Eq. (|3lD, with an epidemic threshold given by 

A - ^^^ - ^ r32l 

In order to compare with the analytical predictions we have carried out large scale simula- 
tions of the SIR model in the WS network with p = 1. In our simulations we consider the 
WS network with parameter m = 3, which corresponds to an average connectivity (k) = 6. 
Simulations were implemented on graphs with number of nodes ranging from A^ = 10^ to 
A^ = 3 X 10^, averaging over at least 10^ different epidemic outbreaks, performed on at 
least 10 different realization of the random network. In Fig. 0, we show the total density of 
removed nodes at the end of the epidemic outbreak as a function of the parameter A. The 
graph exhibits an epidemic threshold at Ac = 0.184(5) that is approached with a roughly 
linear behavior by -Rqo- A linear fit to the form R^o ~ (A — Ac)^ provides an exponent 
f3 = 0.9(1), in reasonable agreement with the analytical finding. This confirms that the SIR 
model in exponentially bounded complex networks has a behavior similar to that obtained 
with the homogeneous mixing hypothesis. Actually, since the connectivity fluctuations are 
very small in the WS graph ((A;^) ~ (k)), as a first approximation we can consider the WS 
model as a homogeneous one in which each node has the same number of links, fc ~ (k). In 
order to provide further evidence to this effective homogeneity, we show in Fig. ^ the time 
evolution of the density of infected nodes for epidemic outbreaks starting only on nodes with 
a given connectivity k. The total epidemic incidence is almost constant for all connectivity 
k, with a slight shift of the peak time of the outbreak. The figure clearly shows that the 
system reacts almost identically to this heterogeneous initial conditions, confirming that the 
homogeneity assumption is correctly depicting the system's behavior. We shall see in the 
next section that this is not the case for SF networks. 



V. POWER-LAW DISTRIBUTED NETWORKS: THE BARABASI-ALBERT 

MODEL 

The Barabasi- Albert (BA) graph was introduced as a model of growing network (such 
as the Internet or the world-wide-web) in which the successively added nodes establish links 
with higher probability pointing to already highly connected nodes [Q. This is a rather 
intuitive phenomenon on the Internet and other social networks, in which new individuals 
tend to develop more easily connections with individuals which are already well-known and 
widely connected. The BA graph is constructed using the following algorithm [§]: We start 
from a small number tjiq of disconnected nodes; every time step a new vertex is added, with 
m links that are connected to an old node i with probability 

m^) = ^, (33) 

where fcj is the connectivity of the i-th node. After iterating this scheme a sufficient number of 
times, we obtain a network composed by A^ nodes with connectivity distribution P{k) ~ k~^ 

9 



and average connectivity (fc) = 2m (in this work we will consider the parameters rriQ = 5 
and m = 3). Despite the well-defined average connectivity, the scale invariant properties 
turn out to play a major role on the physical properties of these networks (for instance, the 
resilience to attack |[I3|JI^ ). 



In the continuous k approximation, that substitutes the discrete variable k for a contin- 
uous variable in the range [m, oo[, the connectivity distribution of the BA model takes the 
form 

2m'^ 
P(k) = — — for k>m. (34) 

k-^ 

For this distribution, the first moment is finite, (k) = 2m, but the second moment diverges 
with the network size, (k"^) ~ logA^. In view of the general result, Eq. ([20|), we observe 
that the epidemic threshold in this particular network tends to zero for large A^. Also, the 
general solutions obtained in Sec. 4 cannot be applied, and we must work out the particular 
solutions of Eqs. (^6|) and ([T5|). 



The equation for R^, with the connectivity distribution (0) is 

i?oo = 1 - 2m2 / k-^e'^^'^^dx = 1 - 2zM x-^e-^dx, (35) 

Jm J z 

where we have defined the new variable z = Xmcf). Performing the integral, we obtain 

R^ = l-e~\l-z)-z^T{Q,z), (36) 



where r(a, z) is the incomplete Gamma function ^\. For small values of z, the function 



r(0, z) can be expanded in the expression [R7 



r(0, z) ^ -(7s + ln(z)) + z + 0{z^), (37) 

where 7^; is the Euler's constant. By inserting this expansion into the expression for R^o, 
we obtain for small values of (p^o 

R^c^2z = \ ik) 0OO. (38) 

On its turn, the equation for 0(t), with the connectivity distribution ( |3^ is 

'^'^^^^ 1 - <p{t) -m r k-^e-^'^^dk. (39) 



dt 
Defining the new variable z = Xm,(f), we can rewrite the previous equation as 

I dz ^ z /""^ _2 -X 



1 - z x-^e-'^dx. (40) 



Xm, dt \m 

In order to study the limit z{(f)) — > we must first integrate by parts the integral in Eq. (^ 
to get 



1 dz 



z f°° 

1-- e~''' - z\a{z)e~''' - z j \a{x)e~''dx. (41) 



Am dt 

10 



We can now take the limit z ^ in the last integral, obtaining 

/ \n{x)e~'^dx = -7^;, (42) 

Jo 

where again 'Je is Euler's constant. Introducing this approximation, and Taylor expanding 



the expression (|T), we obtain 

1 dz 



'-^^-A^-^"^^^ 



Xm dt 
This equation can be integrated, to yield 

1 /. 1 



(43) 



1-^ 



^exp 1 - 7e - ^ - Ae-^™* , (44) 



Xm \ Xm 

where A is an integration constant. The stationary regime for long times is 



e 



I-JE 



■'00 — 



, -e-^/^'", (45) 

Am 

and by inserting this result into the expression for the total epidemic incidence we find 

Roo ~ e-'/^'". (46) 

That is, the function Roo is non-zero for any non-zero value of A, which is in agreement with 
the predicted threshold Ac = 0. This result also recovers the same behavior obtained by 
considering a diverging connectivity variance in the results reported by May and Anderson 
for HIV spreading in heterogeneous population [0,^. 

The numerical simulations performed on the BA network confirm the picture extracted 
from the analytic treatment. We consider the SIR model on BA networks of size ranging from 
A^ = 10^ to A^ = 10^, with m = 3 and thus (k) = 6. As predicted by the analytic calculations. 
Fig. H shows that R^o decays with A as Roo ~ exp(— C/A), where C is a constant. In order to 
rule out the presence of finite size effects hiding an abrupt transition (the so-called smoothing 



out of critical points ^7j), we have inspected the behavior of the stationary persistence for 
network sizes varying over three orders of magnitude. The total absence of scaling of i?oo 
and the perfect agreement for any size with the analytically predicted exponential behavior 
allows us to definitely confirm the absence of any finite epidemic threshold. A closer look at 
Roo is given in Fig. ^. While Fig. ^ reports the average over 10"^ — 10^ epidemic outbreaks. 
Fig. ^ reports an illustration of the behavior of the cumulative probability P{Roo > R) of 
having an outbreak which affects R individuals in a single realization at A = 0.09. The 
figure shows a finite probability of having outbreaks involving a number of individuals of the 
order of the network size A^. The large plateau corresponds to a gap between large events 
and small outbreaks that give rise to a zero density of infected individuals in the A^ — *■ 00 
limit. Accordingly to the predictions, the plateau extends proportionally to A^ for increasing 
network sizes. 

The large heterogeneity of these networks can be pictorially characterized by inspecting 
the epidemic evolution in each class of connectivity k. We know from Eq. (P^), that the 
susceptibles densities Skit) decay much faster in the highly connected classes. In particular, 

11 



we have that 5*^(00) ~ exp(— Afc^oo)- In Fig. ^, we report Sk{oc) as a function of fc on a 
semi-logarithmic scale. The plot shows the expected linear relation in k. Since Rk{oo) = 
1 — Sk{oo), the curves clearly show that the higher is the nodes' connectivity, the higher is 
the relative incidence of the epidemic outbreak. Classes of nodes with few connections have 
a small density of removed individuals (total number of infected individuals), while highly 
connected classes {k ^ 100) are almost totally affected by the infection. A further striking 
evidence of the peculiar behavior of the SF networks is obtained by inspecting epidemic 
outbreaks starting on nodes with different connectivity k. While an analytical solution 
for this case is very troublesome, the numerical investigation presents clear-cut results. In 
Fig. ^ we present the infection incidence profile for epidemic outbreaks started on sites with 
different connectivities k. The population results much weaker (higher number of infected 
individuals) to epidemics starting on highly connected individuals. This weakness points 
out that the best protection of these networks can be achieved by targeted immunization 
programs [^^. 

VI. GENERALIZED SCALE-FREE NETWORKS 

Recently there has been a burst of activity in the modeling of SF complex networks. The 

recipe of Barabasi and Albert [Q has been followed by several variations and generalizations 

38|^I1 and the revamping of previous mathematical works [Q . All these studies propose 



methods to generate SF networks with variable exponent 7. The analytical treatment pre- 
sented in the previous section for the SIR model can be easily generalized to SF networks 
with connectivity distribution with 7 > 0. Let us consider a generalized SF network with a 
normalized connectivity distribution given by 

P{k) = (1 + 7)mi+^A;-2-^, (47) 

where we are approximating the connectivity k as a continuous variable and assuming m 
the minimum connectivity of any node. The average connectivity is thus 

/•oo 1 + ^ 

{k) = / kP{k)dk = -—^m, (48) 

Jm 7 

while the connectivity fluctuations are given by 

.2\ 7 + 1. 



k') =^—m' if 7> 1, (49) 

/ 7 — 1 



'k^) =00 if 7 < 1. (50) 

Thus, according to the general result Eq. (pO|), the epidemic threshold, as a function of 7 is 



1 



(51) 
(52) 

To obtain the explicit expression for (f)oo and i?oo we must solve the Eqs. (|15]) and (|16D 



Ae = 


7m 


it 


7 


>1, 


Ae = 


= 


if 


7 


< 1. 



for the general connectivity distribution (|47|) . While the differential equation ([15D cannot 
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be solved in a closed form for general 7, we can obtain approximation to the steady state 
value at long times, (poo, solving the algebraic equation 

<^oo = l-7^E^^(^K''^°°- (53) 

\'^) k 

As a function of cpoo from Eq. (|T^), the total epidemic incidence R^q takes the form, 

/•OO /"OO 

1 - i?oo = (1 + iW^^ / dkk'^-^e-^'"^^ = (1 + -f)z^+^ / x-^'^e-'^dx (54) 

Jm J z 

= (l + 7)^'+^r(-l-7,z), (55) 



where we have defined z = Afc^oo and T{a, z) is the incomplete Gamma function [0. In the 
limit 2 — > 0, we can perform a Taylor expansion of the incomplete Gamma function, with 
the form p7| 

r(a, z) = r(a) - - + ;2» V / \ , , (56) 

^ ^ ^ ^ a ^1 (a + n)n! ^ ^ 

where r(a) is the standard Gamma function. This expansion has obviously meaning only 
for a 7^ —1, —2, —3, . . . Thus, integer values of 7 must be analyzed in a case by case basis. 
Substituting the last formula into the expression for Rao-, we are led to 

Ro. = r(-7)z^+^ - (1 + 7) E f u ^ ^ ^^^' + ^(^'^')- (^^) 

;^Jn - 7 - l)n! 7 

That is, for any value of 7 > 0, the leading behavior of i?oo is 

7 + 1 

-Roo ~ 2 = A (A;) 0OO, (58) 

7 



which is equivalent to the expression found for regular networks in Eq. (pH]). 

In order to find the infinite time limit value (poo, we must solve the Eq. (^). Substituting 



the form of the generalized connectivity distribution (^Tj) , we have 



/•OO /"OO 

1 - 7m^ / k-^-^'e-^'l'°°^dk = 1 - 7^^ / x-^-^e-^'dx = 1 - 7^^r(-7, z) (59) 



where again z = \k(poo- Inserting in this last expression the Taylor expansion for the 
incomplete Gamma function, we obtain 

^. = .^r(i-oO + 7i:(i^ (60) 

The leading behavior of the r.h.s. of this equation depends on the particular value of 7 
considered. 

(a) < 7 < 1: In this case, we have 

0oo^(Am0oo)^r(l-7), (61) 
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from which we obtain 

(j)^ ~ (Am)^/(^-^)r(l - 7)i/(^-^). (62) 

Combining this result with Eq. (|58|), we obtain 

Roo ~ A^/(i-^). (63) 

In this range of values of 7 we recover the absence of the epidemic threshold, and the 
associated critical behavior, as we have already shown in Sec. 0. 

(b) 1 < 7 < 2: To obtain a nontrivial solution for 0oo, we must keep the first two terms 
in the Taylor expansion (|60|), namely: 



</.„, ~ (Am0oo)''r(l - 7) + ^— Am0oo. (64) 

7-1 



From this equation we find 

1 1/(7-1) 



which yields 



with an epidemic threshold Ac given by Eq. (^TJ). 

(c) 7 > 2: The most relevant terms in the expansion of 0oo are now 



00 


7 TTl / 7 — l\ 

_r(2-7)(Am)^ V 7^^ J J 




i?^^(A-A,)V(^-^), 



(65) 



(66) 



^00-7^-7^^^. (67) 

7 — 1 7 — 2 

The relevant expression for 0oo is 

^„-l^^fA-^. (68) 

7 — 1 A^m y 7m J 

that for the epidemic incidence yields the behavior 

i?,, ~ (A - Ae). (69) 



The threshold Ac is again given by the general expression (|5T|). In other words, we recover 
the usual epidemic framework in networks with connectivity distribution that decays faster 
than k to the fourth power. Obviously, an exponentially bounded network is included in 
this last case. 
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VII. CONCLUSIONS 

The presented results for the SIR model in complex networks confirm the epidemiological 
picture proposed in previous works. The topology of the network has a great influence in the 
overall behavior of epidemic spreading. The connectivity fluctuations of the network play a 
major role by strongly enhancing the infection's incidence. This issue assumes a particular 
relevance in the case of SF networks that exhibit connectivity fluctuations diverging with 
the increasing size A^ of the web. SF networks are therefore very weak in face of infections 
presenting an effective epidemic threshold that is vanishing in the limit N ^ oo. In the case 
of the SIR model in an infinite population this corresponds to the absence of any epidemic 
threshold below which major epidemic outbreaks are impossible. These results strengthen 
the epidemiological framework for complex networks reported for the SIS model ITHM and 



proposed as well for the SIR model [^. The emerging picture is likely going to stimulate 
the re-analysis of several concepts of standard epidemiology such as the "core group" or the 
characteristic number of contacts that appears to be ill-defined in SF networks. 

The high heterogeneity of SF networks finds signatures also in the peculiar susceptibility 
to infections starting on the most connected individuals and the different relative incidence 
within populations of varying connectivity k. It is reasonable to expect that these features 
can point at better protection methods for these networks which appear to have practical 
realization in many technological and biological systems. In this perspective, the introduc- 
tions of many elements of realism and a better knowledge of the networks temporal pattern 
are fundamental ingredients towards a better understanding of the spreading of information 
and epidemics in a wide range of complex interacting systems. 
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FIGURES 
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A -A, 

FIG. 1. Total density of infected individuals i?oo as a function of A — Ac for the SIR model in 
WS networks of size N = 10® . The value of Ac = 0.184(5) is in good agreement with the analytical 
prediction. The full line is a fit to the form Roo -^ (A — Ac)'^ with (3 = 0.9(1). 
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FIG. 2. Total density of infected nodes as a function of time for the SIR model in WS networks, 
starting from initial conditions peaked in nodes of connectivity k (initial infected individuals ran- 
domly distributed only among the nodes of connectivity k). The spreading rate is fixed to A = 0.20. 
The network size is iV = 10®. 



0:5 




FIG. 3. Total density of infected individuals i?oo as a function of 1/A for the SIR model in BA 
networks of size A^ = 10^ . The linear behavior on the semi-logarithmic scale proves the stretched 
exponential behavior predicted by Eq. (p6[). The inset show the time profile of the average density 
of infected individuals at the spreading rate A = 0.9. 
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FIG. 4. Cumulated outbreak epidemic distribution for the SIR model in BA networks. The 
spreading rate is fixed to A = 0.09. 
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FIG. 5. Density Sfc(oo) of susceptible nodes as a function of k for the SIR model in BA networks. 
Epidemics start from random initial conditions (initial infected individuals randomly distributed 
among all nodes). The spreading rate is fixed to A = 0.09. The network size is A^ = 10^. The 
linear-log plot recovers the exponential form predicted in Eq. (0) 
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FIG. 6. Total density of infected nodes as a function of time for the SIR model in BA networks, 
starting from initial conditions concentrated in nodes of connectivity k (initial infected individuals 
randomly distributed among the nodes of connectivity k). The spreading rate is fixed to A = 0.09. 
The network size is A^ = 10^. 
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